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Abstract: The ejection of material during laser ablation gives rise to the 
O I development of discontinuities in the ambient gas. Several of these discon- 

•^ tinuities are observed and characterized, including externally and internally 

pg propagating shock waves, contact surface, and the ionization front. Quali- 

T-H tative experimental observations and analysis of these discontinuities is pre- 

^_^ sented. Results from shadowgraphy enabled determination of an irradiance 

t^ threshold between two different ablation mechanisms, and determination of 

• ^H several stages of plasma plume evolution. Consideration of the refractive 

Q^ index as a dynamic sum of the contributions from gas and electrons led to 

^ separate identification of ionization front from the contact surface. Further- 

more, ionization front was observed to lead the shock wave at the earlier 
stage of the ablation. 
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1 Introduction 



> 

qn^ Laser ablation is the term used for the process of removal of a small amount 

^L^ of material from its surface (usually, of the order of a few tens of ng to a 

CO /ig) by means of a focused laser beam. Results of extensive studies of laser 

■^ ablation can be found in the textbook authored by Bauerle [ ], while a wide 

^r^ range of applications are discussed in the books [2-4] . Due to complexities of 

. . the processes during ablation, the exact description of ablation mechanisms 

.!_H at different irradiations for different materials remains a difficult problem. 

rS This problem is usually attacked by determining threshold points in the 
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transition of the behavior of the ablation as a function of irradiance (or, 
equivalently, fluence). The range between two threshold points would de- 
fine a single dominant mechanism for ablation, which is then modeled by 
analytical or computational methods. A complication arises when ablation 
involves plasma that interacts with the incoming laser beam [^] . 

From the experimental point of view, ablation mechanisms are deduced 
from the results obtained by several techniques, which can be divided into 
two general types: techniques that investigate the material surface after ab- 
lation, and those that investigate the dynamics of the ablation plume. The 
outcome of the first type is the volume of the ablated material or the abla- 
tion rate, which can be measured by confocal microscopy surface profiles [6], 
or the morphology of the ablated area, which can be determined by scanning 
electron microscopy [7]. The second type of experiments includes visualiza- 
tion and spectroscopy techniques, which measure the shape and velocity of 
the ablation plume, density and temperature of the plume, and variation of 
plume composition, all these resolved in space and time. An example of the 
application of these latter techniques is the study of the ablation process for 
lithotripsy [8-11]. The interest in plasma generation and shock propagation 
in laser lithotripsy is that the mechanisms of stone destruction and removal 
are not yet fully understood, in particular the removal of ablated material 
from the bulk of the stone and the coupling of opto-mechanical properties 
for efficient stone destruction. 

Ablation is invariably coupled with the rapid ejection of material, which 
leads to formation of several dynamic zones in the surrounding ambient. 
This paper presents a study of these zones, with emphasis on the observation 
of discontinuities and their evolution. Since the approach is qualitative, 
only a simple analytical model is taken into consideration, supported by 
experimental results. Shadowgraphy is a well suited and one of the most 
accessible techniques for visualization of these discontinuities. In particular, 
the rapid ejection of the plasma plume and the subsequent evolution of 
discontinuities in the ambient gas easily can be observed with relatively 
good spatial and temporal resolution, while densities and temperatures can 
be assessed only qualitatively. 
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Figure 1: Qualitative illustration of the change of refractive index. 

The number of publications on laser ablation that rely on shadowgraphy 
is enormous, but as an example of the use of the technique might serve the 
well known papers by Callies et al. [I'-] and Brietling et al. [i ]. Similarly to 
these references, our paper reports the observation of discontinuities due to 
ablation, namely externally and internally propagating shock waves, contact 
surface, and ionization front. But, unlike most published works, our paper 
draws attention to the observation that, in a certain stage of the evolution 
of the plasma plume, the ionization front leads the shock wave. 



2 Refract ivity of plasma 

For the ablation laser energies used in this study (irradiation range / = 
6 — 130 GW/cm^), material removal always is coupled with the creation of 
the plasma. In other words, the ablated plume consists of free electrons, 
ions, and neutral atoms. Since the process is dynamic, the number densities 
of these constituents change with time, and at later stages (order of ^s), 
even clusters and congregated microparticles are formed [I I]. In addition, 
during the entire ablation process, the plume interacts with the surrounding 
air, causing the change of its density, mainly through the effect of the shock 
waves [15]. 

As it is well known, shadowgraphy registers the displacement of light 
due to gradients in the refractive index, yielding images with zones of vary- 
ing brightness. The inhomogeneities in the medium through which light is 



propagating cause the light to bend, and the angle of deflection e is related 
to the refractive index n(r) by equation [16] 

ude = Vlnn{r)ds, (1) 

where ds is the infinitesimal path that light traverses and u is the unit 
normal vector. Refractive index of plasma is a combination of the refractive 
indices of its constituents, i.e. neutral gas, ions, and electrons [l ]. As it is 
discussed by Hermoch [ ] , in a first approximation of the dispersion theory, 
the refractive indices of a gas consisting only of neutrals and a gas consisting 
only of ions, are considered to be of the same order of magnitude. The overall 
contribution of neutral and ionic gases on the refractive index of the plasma 
is linearly proportional to the concentration of constituents, i.e. density. 
Therefore, their combined refractivity is expressed by the Gladstone-Dale 
relation 

(n - 1)^ cc p, (2) 

which holds for a wide range of densities. 

On the other hand, the refractive index for the electron gas is 

where Up is the plasma frequency, lo = c/A„ is the frequency of the visual- 
ization light, A'^e is the electron number density, and 

N. = ^ (4) 

is the cutoff density (c speed of light, rUe electron mass, eo vacuum permit- 
tivity, e electron charge). Equation (3) shows that when 

Ne > iVc, (5) 

rie is imaginary, and as a consequence, the visualization light is reflected by 
plasma, and only the near field evanescent wave manages to pass through. 



As the electron number density decreases, plasma transmittance increases 
and equation (3) can be approximated by relation 

(n - 1)^ oc -Ne. (6) 

The overall refractive index of the plasma, then, can be expressed as 

(™ - ^)plasma = (« " l)p + (^ " l)e • (7) 

Obviously, this equation is not an exact one, but it can be useful in demon- 
strating the trend of the change of the refractivity of plasma. In this simple 
case, the two-gas model is applied to plasma and influences from density 
gradient and electron number density variation can be determined. 

Substitution of equations (2) and (6) into (7) shows that refractivity of 
plasma is a contribution of two opposing effects on the angle ue, as expressed 
by (1). Local increase of gas density p has the effect that it increases n 
locally, producing a gradient and bending the light in the direction of higher 
density. This is the familiar case of observation of shock waves, which are 
recorded as adjoined narrow dark and bright lines in shadowgraph images. 
The effect of electrons, though, is opposite: the increase of electron number 
density reduces n, thus bending the light away from the zone with high A'^g- 

The effect of the main discontinuities on the refractivity of plasma is 
illustrated in figure 1, where density jump due to shock wave is modeled on 
a simple blast wave, and where a linear decrease of electron number density 
with distance r from the ablation center is assumed (a study of several profile 
functions for electron number density can be found in [19]). The effect on 
the refractive index is the generation of two discontinuities with gradients 
of opposite sign, which can be recorded by shadowgraphy, as demonstrated 
by the following experiment. 



3 Experiment 

The experiment was done by ablating metallic targets (aluminum and cop- 
per) in air using a laser beam with Xa = 1064 nm, r^ = 7 ns and Ea = 
17 — 360 mJ, focused onto a spot with diameter d = 200 ;um. The choice of 
metallic targets for this study is obvious, because they readily provide an 
excess of free electrons. Placement of the sample was chosen such that air 
breakdown, judged by the luminous volume of the plasma in the absence of 
the target, would be a few millimeters behind the position designated for 
target placement. 

The visualization was done by a standard shadowgraph system, consist- 
ing of a laser beam with A„ = 532 nm and r„ = 3 ns used for illumination, 
and an ICCD image sensor (1280 x 1024 pixels, pixel size 6.7 x 6.7 ^um^) for 
image recording. The visualization beam was expanded and collimated by a 
diverging-converging lens combination, and directed parallel to the material 
surface, perpendicular to the ablating laser beam. A single converging lens 
was used for imaging. The expanded beam had a diameter about two times 
larger than the diagonal of the image sensor. 

Visualization presents two challenges: since plasma emits a very intense 
light, a light source that overpowers plasma's self-luminosity is required. On 
the other hand, the camera sensor should not be overexposed. Our approach 
was to use a combination of a laser beam as a powerful light source with 
a series of neutral density filters as adjusters of illumination that reaches 
the image sensor. Optical density of ND filters was ~ 2.7 (transmittance 
~ 0.2%), and they were placed just in front of the camera, so that the 
entire plasma and ambient light is blocked. Using laser as a light source also 
provided the benefit of having extremely short exposures. The visualization 
laser light had a pulse energy of 25 mJ, but since it was expanded to a large 
area, its influence on the plasma was negligible. 

Scattered light from the exit port of the ablation laser was collected by a 
high-speed biased photodetector (rise time 1 ns) and used as the triggering 
signal for the entire system. Triggering was controlled by a pulsed delay 
generator and was applied to the visualization laser and the camera. 




(a) At = 40 ns (b) At = 180 ns 

Figure 2: Discontinuities from ablation of aluminum by 17 mJ laser energy. 







(a) At = 40 ns (b) At = 80 ns (c) At = 140 ns 

Figure 3: Same as in figure 2, but now laser energy is 25 mJ. 



4 Results 

Two examples of visualization results are presented in figures 2 and 3, which 
show the discontinuities due to ablation of aluminum by irradiances of I = 6 
GW/cm^ and / = 10 GW/cm^, respectively. An image frame corresponds to 
the field of view of 2.1 x 2.1 mm^ (2.1 jivn per pixel), and the time when the 
image is taken is noted in the caption. Time At is defined as the time that 
has passed from the moment when the first motion of the material is observed 
to the moment when the pulse of the visualization laser is fired. It should 
be noted that the initial time of the ejection of the material depends slightly 
on the ablation laser energies, which might come about due to different 
mechanisms that lead to ablation: evaporation and normal boiling, or phase 
explosion [20,21]. 



Figures reveal the disadvantage of using a laser beam for visualization: 
numerous interference fringes that arise from the neutral density filters and 
protective windows attached to the image sensor produce images of aver- 
age quality. Nevertheless, discontinuities are clearly observed. Since the 
visualization laser can be triggered to emit only a single pulse, each image 
corresponds to a fresh ablation. Comparison of images obtained under the 
same conditions reveal that the ablation was done with excellent repeata- 
bility. Images for ablation laser energies higher than 25 mJ show the same 
features as those of figure 3. Similar results are obtained for copper as well. 

Figure 2 shows the familiar sight of a shock wave in the form of a blast, 
expanding away from the ablation zone. The visualization of the shock 
shows a leading dark stripe followed by the bright stripe. No plasma plume 
is observed in this series of figures. Figure 3, on the other hand, is much 
more complex. Initially, as seen in figure 3a, visualization of ablation is 
characterized only by an opaque hemispherical plasma plume. At a later 
stage, (figure 3b) some morphology of the plasma plume can be observed. 
The visualization of ablation is led by a white stripe of light, followed by a 
dark one. The transmittance of the plasma plume inside these stripes varies, 
with a semitransparent and an opaque zone clearly identifiable. At a much 
later stage, as shown in figure 3c, the leading stripe is dark, immediately 
followed by a bright stripe. This is the familiar effect of a shock wave 
in shadowgraphy. Following is a band of bright stripe, and the interior 
of the plasma plume now shows even more complex morphology: plasma 
transmittance has increased, and plasma plume and internally propagating 
shock waves have become visible. At much later times, the visualization 
becomes qualitatively identical to that of figure 2, with only the shock wave 
propagating away from the ablation zone being visible, until the shock is 
dissipated into a sound wave. 

5 Discussion 

The dramatic difference between the discontinuities visualized in figures 2 
and 3 is an indication that the ablation mechanisms were different, but 



shadowgraphy is very limited in offering an explanation. Still, based on 
discussions from Section 2, several observations can be made. 

According to equation (4) , the cutoff electron number density for our ex- 
periment is Nc ~ 10^'' cm~^. For the ablation shown in figure 2, this number 
is not achieved, because plasma is transparent throughout the process and 
only the shock wave is observed. 

The next level of irradiance tested, achieves the cutoff value close to 
the target and for initial moments of plume ejection, as it can be inferred 
from figures 3a and 3b, where plasma is opaque. Further evidence for the 
large electron number density is provided by figure 3b. Based on equation 
7, the refractivity of plasma is dominated by the contribution of electrons, 
because the light is being deflected away from the target. The effect of 
electrons is observable even at later stages (as in figure 3c), although now 
shock waves have expanded and precede. Since the lower limit of electron 
number density that would effect the refractivity of plasma is of the order of 
10^^ cm~^, it can be concluded that N^ is greater than this value. Order-of- 
magnitude comparison to spectroscopic measurements further support this 
conclusion [22]. 

Ablation studies by visualization make use of the blast wave theory [23], 
where the pressure and temperatures involved in the ablation process are 
calculated from the measurements of shock radius [ - ] . The evolution of the 
shock radius is described by equation 

/ P X 1/(2+0 

r(^t) = c(^^^] t"^^+^\ (8) 

V po / 

where C is a constant (usually C ~ 1), Eesw is the energy content of the 
shock wave, and po is the density of the surrounding atmosphere. ^ is a 
parameter that accounts for the dimensionality of the propagation, and it 
is equal to 2 for spherical propagation of the shock wave. Spherical symme- 
try of the propagation of the shock wave is an evidence that plasma being 
initiated by the breakdown of the solid, rather than the breakdown of the 
surrounding air, because, in the latter case, plasma would appear elongated 
along the Rayleigh range, exhibiting a cylindrical symmetry. 
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Figure 4: External shock front velocity as a function of time. 



Results for shock velocities as a function of time are shown in figure 4. 
Points show the velocity calculated as the difference of the radius of the 
shock between two successive images, divided by the time interval between 
those images. The curve is the derivative of r(t) from equation (8). Since the 
energy content of the shock wave does not change much with the increase 
of irradiance, a single velocity curve is shown for illustrative purposes. 

Applying deduced shock velocity to strong shock relations, we can easily 
calculate the pressure and temperature behind the shock wave [2-"']. The 
relation between pressure and shock velocity is given by 



Po 



drV 



(9) 



-f + l\dtj ' 
where 7 is specific heat ratio of the surrounding ambient (for air, 7 = 7/5), 
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Figure 5: Ablation discontinuities and zones. Ablation laser energy used 92 
mJ. At = 180 ns. 



while the relationship for temperature is 

7-1 P 



T = T< 



7 + l^o' 



(10) 



with To and Pq being room temperature and pressure, respectively. Calcu- 
lations show that within the first 300 ns, pressure and temperature of the 
shocked gas just behind the external shock wave decrease from 300 MPa and 
100 000 K down to 0.5 MPa and 2 500 K. Since the shock Mach number 
throughout this time interval is larger then 10, it can be stated that the 
shock wave maintains the plasma by partially ionizing the ambient air. 

In summary, different discontinuities of the ablated plasma plume can be 
defined as in figure 5: ESW - externally propagating shock wave, IF - ion- 
ization front, CF - contact front, ISW - internally propagating shock wave. 
ESW divides the laboratory air and the shocked air, which is then followed 
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by the zone dominated by an excess of electrons - the ionized air. Since ini- 
tiahy the ionization front was observed to be ahead of the shock wave, but 
later passed by the shock, the ionized air is initiated by the electrons from 
the material, but later, the strong shock wave is able to maintain a level 
of ionization. Contact front is the surface that designates the beginning of 
the ablated material - plasma plume. Its constituents are cooled ions and 
neutrals from the target that collide with the air molecules. Since the target 
is continually ejecting hot material, this leads to the creation of the internal 
shock wave, which initially moves away, but later moves towards the target. 
Based on figure 3, a rough estimate of time stages for visualization can be 
set as follows: during the first stage (~ 50 — 60 ns) the number density 
of electrons is larger than the cutoff value, therefore the plasma is opaque, 
because it completely reflects the visualization beam. During the second 
stage (~ 60 — 80 ns), the ionization front leads the shock wave, and plasma 
is semi transparent. During the third stage (« 80 — 160 ns), external shock 
wave overtakes the ionization front, internal shock wave dissipates, and the 
entire ablation process can be modelled as a blast wave. Finally, at later 
stages, only the external shock wave is visualized. 

6 Conclusions 

In this paper, the observation and interpretation of discontinuities following 
laser ablation were discussed based on the change of plasma refractivity. 
It was shown that shadowgraphy provides a simple, but useful method for 
identifying laser ablation mechanisms and for plasma diagnostics, although 
only order-of-magnitude quantitative results can be obtained. Specifically, 
a transition in the ablation mechanisms was observed between irradiances 
6 and 10 GW/cm^. The approximate electron number density was judged 
based on the transmittance of the plasma. 

A more useful result from shadowgraphy was the ability to identify sev- 
eral discontinuities in the ambient surrounding the ablation spot, and mea- 
sure their evolution with good spatial and temporal resolution. This allowed 
for characterization of the ablation plume and definition of zones of different 
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constituents, based on their effect on the refractive index. Special attention 
was drawn to the observation of ionization front leading the externally prop- 
agating shock wave, which indicated that initially plasma was generated by 
ejected electrons from the target, and only at later stages the shock wave 
was able to induce air ionization. 
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